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SUMMARY: In this article we examine the possibility to extend leap second 
extrapolation for a near future based on some periodic terms in the Earth's rota- 
tion changes. The lERS data, covering the interval from 1962.15 to 2006.95, are 
analyzed. The difference AT is extrapolated till to 2035 and compared with the 
lERS extrapolated values to the 2012. It can be seen that for the interval from 2006 
to 2024 only 1 leap seconds (negative) will be operated. 
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1. The leap second 



The recommendations of the lAU (Interna- 
tional Astronomical Union) were formalized by res- 
olutions of their Commissions that the name UTC 
(abbreviation is compromise between English Coor- 
dinated Universal Time and French Temps Universel 
Coordonne) was retained. UTC was recommended 
as the basis of standard time in all countries, the time 
in common (civil) use. The limit of [UTl — UTC] 
{UTl, Universal time) was set at -1-0.950 s, as this 
is the maximum difference that can be accommo- 
dated by the code format. The maximum devia- 
tion of UTl from [UTC + DUTl] (time correction 
DUTl = UTl - UTC) was set at -HO. 100 s. In 
1974, the CCIR (Consultative Committee on Inter- 
national Radio, En. or Comite consultatif interna- 
tional pour la radio, Fr.), increased the tolerance for 
[UTl - UTC] from 0.7 s to 0.9 s. The present UTC 
system is defined by ITU-R (International Telecom- 
munication Union - Radio) (formerly CCIR) Recom- 
mendation ITU-R TF.460-5: 

"UTC is the time scale maintained by the 



BIPM (Bureau International des Poids et Mesures, 
Fr. or International Bureau of Weights and Mea- 
sures, En.), with assistance from the lERS (Inter- 
national Earth Rotation and Reference Systems Ser- 
vice) , which forms the basis of a coordinated dissem- 
ination of standard frequencies and time signals. It 
corresponds exactly in rate with TAI (International 
Atomic Time, En., Temps Atomique International, 
Fr.) but differs from it by an integral number of sec- 
onds. The UTC scale is adjusted by the insertion 
or deletion of seconds (positive or negative leap sec- 
onds) to ensure approximate agreement with UTl." 

Because the Earth's rotation is gradually slow- 
ing down, and in addition it has both random and 
periodic fluctuations, it is not a uniform measure of 
time. The time difference AT = [ET - UTl] = 
[TT ~ UTl] represents the difference between the 
uniform scale of Ephemeris Time (ET) or Terres- 
trial Time (TT) and the variable scale of Univer- 
sal Time (UTl). Before 1955, the values are given 
by AT = [ET - UTl] based on observations of 
the Moon. After 1955, values are given by AT = 
[TT-UTl] = [rA/ + 32M84-C/Tl] from measure- 
ments by atomic clocks as published by the BIH (Bu- 
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reau International de I'Hcurc, Fr. or International 
Time Bureau, En.) and the BIPM. 

According to Stephenson and Morrison 
(1995), over the past 2700 years AT can be repre- 
sented by a parabola of approximately the form 

AT = (3l7cy^)(T - 1820)^/(100)^ - 20^ 

where AT is expressed in seconds and T is the year. 
The minimum is at the year 1820 and passes through 
at the year 1900. Actual values of AT based on 
astronomical data may differ somewhat from this 
smoothed fit. For example, the value of AT is 
32M84 at 1958.0, the origin of TAI. However, no 
single parabola can satisfactorily represent all the 
observational data. 



The derivative of AT is 



ALOD = (0.0017 sd-V"^) 



(T - 1820) 

Too ' 



which represent the change of the length of day 
(LOD) in SI seconds. Different studies implies differ- 
ent values. The actual value of the LOD will depart 
from some long-term trend with short-term fluctua- 
tions (periodicity) between ±3 ms on time scale of 
decades. 

Similarly, the insertion of leap seconds is due 
to the fact that the present length of the mean solar 
day is about 2.5 ms longer than a day of precisely 
86400 SI seconds, as a consequence of the long-term 
trend, so that the Earth's rotation runs slow with 
respect to atomic time. 



2. The data and model 



UTC is kept within 0.9 s of UTl by the occa- 
sional insertion of a leap second adjustment. When 
the present UTC system was established in 1972, the 
time difference AT = [TT-f7Tl] = [TAI+Z2\1M- 
UTl] was equal to 42.23 s. Thus the difference be- 
tween TAI and UTl in 1972 was approximately 10 s. 
To maintain continuity with UTl, UTC was initially 
set behind TAI by this amount. As of January 1st, 
2006 the 23 positive leap seconds have been added. 
Thus UTC is presently behind TAI by 33 s. Figure 1 
illustrates the relationships between TAI, UTC and 
UTl (lERS data). 

A least-squares fit of the difference [UTC — 

TAI] since 1972, shown in Figure 1, implies a nearly 
linear increase with a slope of (2.10 + 0.05) ms per 
day. This value represents the average excess in the 
length of day during the past three decades and is in 
approximate agreement with the value computed on 
the basis of the long-term trend. 

Recent global weather conditions have con- 
tributed to a short-term change in the length of day. 
Decade fluctuations due to the interaction between 



the Earth's core and mantle and global ocean circu- 
lation may also contribute. The model characterized 
by triaxial Earth's structure, its fluid core, visco- 
elastic mantle and equilibrium ocean was proposed 
by the Vondrak (1987). 




1990 2000 

year 



Figure 1. UTC - TAI as extrapolated by lERS 
at 2001 and 2002 compared with observed values of 
UTC - TAI till 2006. 

As a contrary Yoshido and Hamano (1993) 
have proposed that one of the main causes of the 
secular variation of the geomagnetic field is length- 
of-day (LOD) variation, namely, the variation of 
the rotation velocity of the mantle. They devel- 
oped an analytical model, in which fluctuations of 
the rotation velocity of the mantle induce flow in the 
outer core through topographic coupling at the core- 
mantle boundary (CMB). The flow in turn bends the 
toroidal field to produce a poloidal field. 

At periods longer than few years, extending to 
many tens of years, the so-called decadal variations, 
the sources of excitation for both LOD and PM are 
more enigmatic. The difficulty is that at these peri- 
ods other effects may be important, including visco- 
elastic behavior such as post-glacial rebound, and 
exchange of angular momentum with the fluid core. 
In particular, it has become common to invoke the 
core as the major cause of decadal LOD changes. 
Although some climatic forcing of long period LOD 
has been recognized (Salstein and Rosen, 1986; Eu- 
banks, 1993), it is uncertain at what time scale air 
and water become less important than the core. Un- 
fortunately, the role of the core remains largely un- 
quantified because it is too remote to be easily ob- 
served. A further difficulty in assessing the air/water 
role at long periods is that the torques required to 
cause decadal LOD variations are utterly insignifi- 
cant when compared with those applied by the at- 
mosphere at shorter time scales (Hide and Dickey, 
1991). This means that atmospheric/oceanic torques 
of geodetic significance are of second-order impor- 
tance in general circulation studies. Quantification 
of momentum budgets among Earth, air, and wa- 
ter reservoirs is thus lacking at long periods. The 
requirements for progress in this field coincide com- 
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pletely with the central problems of global climate 
change. 

The five Earth Orientation Parameters (EOP) 
- two components of polar motion a;,y, two compo- 
nents of celestial pole offsets Aip, Ae, and univer- 
sal time UTl (that is nothing else but the angle of 
Earth's rotation around its spin axis) - are now ana- 
lyzed and routinely derived from the observations at 
several world's centers, combined and regularly pub- 
lished in lERS bulletins. The most recent system of 
constants and algorithms (McCarthy, 1996, Vondrak, 
1998) are used . The past solutions based on optical 
astrometry were merged with the combined solutions 
from the modern techniques (Vondrak, 2001). 

The existing ERF (Earth Orientation Param- 
eters) series have been analyzed by many scientists. 
The most extensive reviews, with historical meaning, 
were given in well-known monographs by Munk & 
MacDonald (1960) and Lambeck (1980). Most recent 
of these analyses is that by Zharkov et al. (1996). A 
very detailed review, mostly concentrated on modern 
space techniques and discussion of short-periodic ef- 
fects, was published by Eubanks (1993). Moreover, 
in our previous paper (Segan et al., 2003) we have 
introduced the relax period as good explanation of 
the residuals in the {UTl - UTC)bli data. 

The motivation for the leap second, therefore, 
is due to the fact that the second as presently defined 
is "too short" to keep in step with the Earth. How- 
ever, had the second been defined to be exactly 
equal to a mean solar second at the origin of TAI 
in 1958, the discrepancy v^rould not have been re- 
moved; the agreement between the SI second and 
the mean solar second would have only been tempo- 
rary and their difference would simply have become 
gradually more apparent over this century. 

Continuing use of a non-uniform time scale in- 
cluding leap seconds in the face of these considera- 
tions could lead to the necessity to proliferate an ef- 
fective method for extrapolation of the future values 
of AT. 

To try that, according to our analysis of the 
lERS AT data from 1962.15 to 2006.95, existence of 
some periodic terms is discovered. By using the har- 
monic analysis 17 components plus a linear term of 
equation (1) are determined. 





+ Aj sin 


r27r \ 













(1) 

In equation (1) A/AT^ is a modified AT^. In order to 
obtain the ATi values one needs to perform a trans- 
lation by a constant value of 32.184 + 10 seconds: 

AT, = 32.184 + 10 - MAT,. 



Table 1. The coefEcients for 17 components ob- 
tained from the harmonic analysis. 
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Figure 2. The dashed line (blue) represents lERS 
data and the solid one (red ) represents our approxi- 
mation (prediction) between 1962.15 and 2006.95. 

The coefficients of the linear function (in equa- 
tion (1)) are: free term Ci = 3.408457 and slope 
coefficient C2 — 0.357042. The coefficients of the 
harmonic components are given in Table 1 . cr| is the 
upper limit for the coefficients of the harmonic com- 
ponents (it is the same for both coefficients, A^ and 
Aj) and it is given with two significant digits. The 
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periods Pj are given in years. Because of the secular 
and decade variations all terms corresponding to in- 
significant periods and amplitudes (smaller than few 
milliseconds) are on the noise level so that they are 
not determined. 

In Fig. 2 the dashed line (blue) represents 
the lERS data from 1962.15 to 2006.95, whereas the 
solid one (red) represents our approximation for the 
period between 1962.0 and 2007.0. One can notice 
an excellent agreement. 

According to our approximation till 2024 there 
is practically no need to introduce leap seconds be- 
cause an accumulation of 0.9 seconds is reached as 
late as at 2014.0, i. e. an accumulation of one second 
in 2015. The local maximum in our approximation 
occurs at 2016.4 and its value is 65.944 which exceeds 
the difference of one second by 0.1 seconds only (Fig. 
3). In view of the further negative trend of our ap- 
proximation the next leap second (negative) should 
be introduced in 2024 if the difference of 0.1 seconds 
at the moment of approximation maximum were ne- 
glected. 
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Figure 3. The approximation from 1990 to 2035 
and lERS data till 2006.95. The indicated points 
are: the moment when the leap second was intro- 
duced for the last time (2006), the next introducing 
(2009), the point at which according to our approx- 
imation the difference of 0.9 seconds between UTl 
and UTC is reached (2014.0), the local maximum of 
our approximation (2016.4), the point for which AT 
is approximately the same as in 2006 when the leap 
second was introduced for the last time (2021.95) 
and the point when, according to our approximation, 
the next leap second (negative) should be introduced 
(2024.0). Note that this means to introduce only one 
negative leap second within 18 years (2006-2024). 
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Figure 4. The historical observations covering the 
period 1620-2006.95 are represented by the dashed 
line (blue) and our approximation concerning the in- 
terval 1620-2035 by the solid line (red). 



4. CONCLUSIONS 



As the day at present is actually closer to 
86400 SI seconds the leap seconds have not been 
required regularly. However, this condition cannot 
persist and the long or mid-term trend will be even- 
tually restored. The asymptotic behaviour of the AT 
in the Stephenson-Morrison approximation (1995) is 
not natural. 

As can be seen from the analysis (Fig. 4) ap- 
plied within a relatively short time interval, 1962.15- 
2006.95, the extrapolation of the AT value agrees 
well with both the historical values (1620-1962) and 
the real measurements from 1962.15 to 2006.95. The 
discrepancy in the amplitude is higher for the for- 
mer period, but the agreement for the phase is very 
good which indicates that the approximation could 
be improved if a bigger set of measurements were 
available. 

Our conclusion, that the first (additional) leap 
second should be introduced in 2014 only, opposes 
to all the extrapolations proposed till nowadays by 
both international institutions and individual ex- 
perts. The majority of them predicts 2008 as the 
year when the additional correction should be intro- 
duced. 

The values attained by AT are of such order 
that all the physical factors unambiguously recog- 
nized till nowadays cannot cause such a phenomenon 
and in this sense the present analysis means a strictly 
mathematical approximation only. We believe that 
the years to come will show the correctness of the 
obtained results. 
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UDK ... 



HcnHTana je MoryliHOCT /la cc npornosnpa 
6poj npecTynHHX ccKynAH y 6jiHCKoj Gy/iyRnocTH 
na OCHOBH no3HaBaH,a hckojihko xapMOHHjcKHX 

KOMnOHCHTH y BcMJLIIHoj pOTaLI,HjH. AnajIHSH- 

paHH cy noAau;H lERS-a kojh oGyxBaiajy HHiep- 
eaji OA 1962,15 ao 2006,95. PasjiHKa AT je eK- 



cTpanojiHpana 2035. ro/jHHC h ynopeljena ca 
npe/iBiil5aH)HMa lERS-a ao 2012. ro/jHHC. MojKe 
ce BH/jcTH je y HHiepBajiy oa 2006. ao 2024. 
rojjHHC noTpe6HO yBecTH caMO je/iny (neraTHBHy) 
npecTynny ceKyuAy- 
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